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I present a brief introduction to the lattice formulation of quantum field theory, and 
discuss the use of lattice simulations for studies in particle physics phenomenology. 
The computation of fs, the decay constant of the B-meson, is used as a case 
study. I also explain the appearance and cancellation of "renormalons" in the 
evaluation of power corrections (higher-twist corrections) in hard scattering and 
decay processes. 



Introduction 

The lattice formulation of quantum field theories, together with large scale 
numerical simulations, is becoming a very powerful non-perturbative tool in 
many areas of particle physics, and in the evaluation of long-distance strong- 
interaction effects in physical processes, in particular. In these lectures I will 
present a brief introduction to lattice computations in QCD (lecture 1) and 
will then discuss, as a case study, the applications to the evaluation of the 
decay constant of the B- meson, and other physical quantities, in the Heavy 
Quark Effective Theory (HQET). Finally in the third lecture 1 will review the 
question of "renormalons" , which arises when one evaluates power corrections 
(i.e. higher-twist corrections) to hard scattering and decay processes. 

The principal objective of these lectures is to explain the theoretical frame- 
work needed for lattice simulations in particle physics, and to give some idea 
of what is possible now and what the outstanding problems are. In the first 
two lectures I will illustrate the discussion with numerical results obtained by 
the UKQCD collaboration, of which 1 am a member. It is not my intention, 
nor would it be appropriate at this school, to present a critical review of the 
many interesting results which have appeared in the last few years from many 
different groups. The proceedings of the annual symposia on Lattice Field The- 
ory contain detailed reviews and status reports on all aspects of the subject, 
and provide very helpful starting points for literature searches. A selection 
of review talks, which are suitable as further reading to the material of these 
lectures, can be found in refs.tlH^. 
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1 Lecture 1 — Lattice Computations of Masses cind Matrix Ele- 
ments 

In this first lecture I will briefly review the ingredients used in evaluating 
hadronic masses and matrix elements in QCD. Simulations using the Heavy 
Quark Effective Theory will be described in the remaining two lectures. 

The quantity which is evaluated directly in lattice computations is the vac- 
uum expectation value of a multi- local operator 0(a;i,a;2, • • • ,a;„), composed 
of quark and gluon fields: 

{0\O{xi,X2,---,Xn)\0) = ^ J[dA^][di;][d^]e'^0{xi,X2,---,Xn) , (1) 
where Z is the partition function 

Z = J[dA^][dilj][chp]e'^ (2) 

and S is the QCD action. In these introductory paragraphs all the formulae 
are written in Minkowski space, although, of course, the lattice simulations 
use a Euclidean formulation of QCD. For the computations described in these 
lectures, n will be either equal to 2 or 3, and we now describe the significance 
of the two- and three-point correlation functions respectively. 



Two-Point Correlation Functions: Consider the bilocal operator O: 

0{XuX2) = J{xi)JHx2), (3) 

where J and are interpolating operators which can destroy or create the 
hadron H, whose mass will be denoted by tuh- Moreover, we will assume 
that H is the lightest hadron which can be created by J^. We now define the 
two-point correlation function: 

C2 {t) = j d^x e'P"'-^(0 1 J{x, t) Jt (0, 0) I 0) , (4) 

where we take t>0. Inserting a complete set of states, {n}, between the two 
operators we obtain 

C2{t) = ^ f d^xe'P-^{0\J{x,t)\n){n\j\6,0)\0} (5) 

n ■' 

= -L|(o|J(0,0)|/f(p))|2e-'^* + ... , (6) 
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where the elhpses represent contributions from states which are heavier than 
iJ, and we have used translational invariance to displace the argument of J 
from (x, t) to the origin. The four- momentum p has components {p, E), where 
E'^ ^ + mjj. In Euchdean space the exponential factor in eq.(^ is e~^*, 
and the contributions of the heavier states have similar exponential factors. 
Thus for sufficiently large times t, the contribution from the heavier states is 
suppressed, and the ground state is isolated. By fitting the time behaviour of 
the computed correlation function, one obtains E, and hence the mass of the 
hadron H, and also the matrix element |(0 | J(0,0) \H{p))\. For example if H 
is the pion, and J the axial current, then from this matrix element we can 
obtain the value of the pion's decay constant: 

\{0\A^iO)\nip))\^Up^. (7) 

Three- Point Correlation Functions: It will also be useful for us to con- 
sider three-point correlation fmictions: 

C3{t^,ty) = I d\ d?y e^P-^e'«-^"(0 | J2(x, t^) 0{y, ty) 4(0, 0) | 0) , (8) 

where, Ji and J2 are the interpolating operators for hadrons Hi and H2 respec- 
tively, O is a local operator, and we have assumed that > ty > 0. Inserting 
complete sets of states between the operators in eq. we obtain 

-iEity „-iE2{t^-ty) 
C3{t.,ty) = — (0 I MO, 0) I H2ip, E2)) X 

{H2{p, E2) I 0(0, 0) I ffi(p + g. El)) {Hi{p + q, El) I 4(0, 0) I 0) • • • (9) 

where Ei = ^mf + {p + q)^, E2 = \/w^^pP- and the ellipses represent the 
contributions from heavier states. In Euclidean space the exponential factors 
become exp(— E'l^j^) and exp(— £'2(^2; — ty))-, so that for large time separations, 
ty and tx — ty, the contributions from the lightest states dominate. All the 
elements on the right-hand side of eq.(|^) can be determined from two-point 
correlation functions, with the exception of the matrix element {H2\0\Hi). 
Thus by computing two- and three-point correlation functions the matrix ele- 
ment {H2\0\Hi) can be determined. 

The computation of three-point correlation functions is useful in studying 
weak decays of hadrons (e.g. if Hi is a D-meson, H2 & K meson and O the 
vector current c'^^s, then from this correlation function we obtain the form 
factors relevant for semileptonic D K decays) and in determining properties 
of hadronic structure (e.g. if Hi and H2 both represent the proton and O is 
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the electromagnetic current, then from C3 we obtain the electromagnetic form- 
factors of the proton). 



The computation of correlation functions will be outlined in subsec. 1.2 
but I start with a brief description of the formulation of QCD on a discrete 
space-time lattice. 



1.1 Elements of Lattice QCD 

In order to perform lattice computations it is necessary to formulate QCD 
in discrete space-time and in Euclidean space ^. The quark fields '^'i^) ^-re 
field variables defined on the sites, {a;}, of the lattice. The gluon fields are 
conveniently introduced in terms of "link" variables C/^ (x) , defined on the link 
between the point x and x + a/i, where fi is the unit lattice vector in the /i 
direction. The gluon field A^, is then given by 



(10) 



where a is the lattice spacing and go (a) is the bare coupling constant defined 
with the lattice action being used. Under a local gauge transformation g{x), 
the fields transform as: 



and 



U^{x)^g{x)U^{x)g\x + afL) . (11) 



Uj^{x—ajl) can be thought of as the link variable from x to x—afi. From eg. (pi]) 
it can be seen that any closed loop of link variables will be gauge invariant, and 
this is exploited in constructing the action and the operators which represent 
observables. For example, the Wilson action for the gauge fields is defined in 
terms of the "plaquette" variable, which is the product of link variables around 
an elementary square of the lattice: 



7'^,,(x) = Tr [U^{x)U,,{x + afi)Ul{x + au)Ul{x)] 
The Wilson action is defined by 

1. 



1- -ReVf,,^{x) 



(12) 



(13) 



where f3 = 6/(7q. Expanding the right hand side of eq.(p^ in powers of the 
lattice spacing, one can verify that, as required, it gives / d^XjF^^, up to 
"irrelevant" terms of O(a^). 



^For excellent pcdagological introductions to the lattice formulation of QCD see refsL 
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When trying to discretise the fermion action, one encounters the famous 
problem of "Fermion Doubhng" . For example consider the free quark propa- 
gator, corresponding to the Dirac action ?/)(^ + m)'0- Defing the derivative of 
a function as the difference of its values at neighbouring sites, one obtains the 
following lattice propagator in momentum space 

'S'(g')freo = m + - ^ 7^ sin(ag^) . (14) 

Consider now q satisfying the continuum mass-shell condition + m? = 0. 
The free-propagator then has a pole at this choice of q as expected (modified 
by terms which vanish as a — > 0). Now there are also poles at the values 
of momenta with any component replaced by -k /a — g^, leading to sixteen 
states instead of one. Such doubling of the fermion degrees of freedom is a 
general feature of lattice theories B. Wilson's solution to this problem was to 
add an "irrelevant" operator to the action, in such a way that the unphysical 
doublers have infinite masses as a ^ 0. For the free theory the Wilson term 
takes the form: 

-yV^aV = -m4^{l^[^(a;)V'(a: + a/i) 

+ Tp{x)'>p{x - afi) - 2lp{x)'>p{x)]^ , (15) 

where r is an arbitrary parameter. Physical quantities must, of course, be 
independent of r. With the Wilson term, the free propagator becomes 

S{q)free = m + - ^jf,sm{aqf,) + - ^{1 - cos{aqf,)) . (16) 

from which the shift of the mass of the doublers can be deduced. The Wilson 
term breaks chiral symmetry, evenior m — 0, but the symmetry is restored 
at a special value of the bare massEI (see below). Including now interactions 
between quarks and gluons, the Wilson action becomes: 

Sw = ■S'giuon + «^ { ~ - [iPix)ir - j^)Uf,{x)ip{x + afi) + 

X fl 

^{x + afi)ir + j^)Ul{x)^ix)] -f V^(x)V'(x)} , (17) 

where the parameter k = l/(2m -I- 8r) contains the dependence on the mass 
of the quark. The fields have been rescaled by a factor of 2k, and of course, 
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a new term must be included for every quark flavour (each with a different k, 
parameter). 

I end this subsection with some brief comments: 

i) The only parameters of the theory are those of QCD itself, i.e. the 
coupling constant, and the quark masses (parametrised in terms of k). 

ii) As mentioned above, physicial quantities must be independent of r, and 
we set r = 1 in most of the following discussion. 

iii) There is a value of k {— Kc) for which chiral symmetry is restored, (kc = 
l/8 + 0(a,))l. 

iv) Expanding the right hand side of eq. (|l^) in powers of the lattice spacing 
we find 

Sw = S^q°S'd + 0{a) , (18) 

where the superscript "cont" implies that this is the continuum QCD 
action. Thus, for a finite lattice spacing a, we expect there to be "dis- 
cretisation" errors of 0(a) (i.e. of O {pa), O {ma) or ofO(AQCDa), where 
p is the momentum of one of the hadrons). Much work is currently be- 
ing done in trying to reduce thjese lattice artefacts, by using "improved" 
lattice actions and operators 



1.2 Evaluation of Correlation Functions 

In this subsection I will briefly describe the steps used in the evaluation of 
correlation functions, without explaining the algorithms which are used in the 
computations. 



Pure Gauge Theories: In the pure gauge theory, and in Euclidean space 
the functional integral in eq.(|l|) becomes: 

{0\O{xi,X2,---,Xn)\0) ^ ^ J[dU]e-^''^^^0{xi,X2,---,Xn) , (19) 
where Z is given by 

Z = J [d[/]e-^(^) (20) 

and S is the gluon action (|l^). Examples of interesting studies in pure gauge 
theories include the evaluation of the spectrum of glueball states and the de- 
temrination of the potential between a static quark and an aiiti-quark (for a 
recent review, and references to the original literature, see ref. til). 
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On a finite lattice with V points, the right hand side of eq.(|l^) is a well- 
defined multi-dimensional integral (specifically, for an SU (n) gauge theory it 
is a AV{n? — 1) dimensional integral). Monte Carlo algorithms are an efficient 
way to evaluate these multi-dimensional integrals. With these algorithms gluon 
configurations are generated on all the links of the lattice, with a probability 
density 1/Zexp(-S'([/)). Thus 

We 

^0)^^Y.O{Uc) (21) 

^ c— 1 

where 0{Uc) is the value of O for the field configuration {Uc}, and the sum is 
over Nc statistically independent gluon configurations. In general we do not 
very know much about the statistical distribution of the results, but expect 
that the uncertainty will typically decrease like 1 / y/Wc. The "statistical" error 
for a quantity in a given computation is estimated from the variation of the 
result as gluon configurations are subtracted and added. 

Lattice QCD We now write the Wilson action ( [TtI ) in the following short- 
hand notation: 

=V^A(C/)^ + ^gluon (22) 

wher S'giuon is the gluon action and A is a matrix in position, spin and colour 
spaces. The functional integral over the quark fields is straightforward to 
evaluate formally: 

y[dV]W]e^«' = det[A(f/)] e-^s'"°" (23) 

J [d^][d^P]4>,{^j e^'^' = A-.i(f/)det[A([/)]e~^— (24) 

etc. The labels i and j represent the lattice coordinates, as well as spin and 
colour indices. A~^{U) is the quark propagator in the background field {U}. 
Thus for the integration over the gluon fields, the action should now be modified 
to S'giuon — In det[A([/)]. The second term, which represents the effects of closed 
quark loops, is non-local, and its inclusion in the Monte-Carlo algorithms is 
very expensive in terms of computing resources. It is the focus for much 
of the development work in improved algorithms for lattice simulations. At 
present, most large scale numerical simulations are performed in the "quenched 
approximation" , in which the determinant is not included, i.e. the effects of 
closed quark loops are neglected. Even where they are explicitly included, it 
is generally the case that the sea quarks are heavy (with masses about those 
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of the strange quark) , so that it is hkely that the effects of hght quark loops 
are significantly underestimated in these studies. 

We do not really know how large the errors introduced by the quenched 
approximation are, and expect that these will depend on the physical quantities 
being computed, as well as on the quantities being used to fix the lattice 
spacing (and quark masses). Over the last ten years and more, a wide variety 
of quantities have been computed in quenched simulations and compared to 
the experimental values. In almost all cases there is good agreement between 
the computed and measured values (of the order or better than 25% or so), 
giving us confidence in the predictions for unknown physical quantities, such 
as the decay constant Jb discussed in the second lecture. The errors due to 
quenching are likely to be larger when one compares physical quantities which 
depend on different scales (e.g. if one fixes the lattice spacing in some process 
which occurs at a scale of O(Aqcd) and uses it in the prediction for heavy 
quarkonia). Of course it is only when we are able to include light quark loops 
in a reliable way that we will be able to claim that lattice QCD is a truly 
quantitative non-perturbative technique. 

For the rest of these lectures I will discuss results obtained with the 
quenched approximation. 

1.3 The Pion Propagator 

I conclude this lecture with a specific example, that of the evaluation of 
the pion's decay constant f^^. The evaluation of other matrix elements and 
hadronic masses follows similar steps. Consider now the correlation function 
of two axial currents: 



where x — {x, t) and 1 have denoted the quark propagator (with the depen- 
dence on the background field implicit) by S{x,Q). This correlation function 
is represented by the diagram in fig. 0. 

The symmetry property S'(0, x) = 7^ S'^{x, 0)7^ is particularly useful, as it 
implies that in order to evaluate the correlation function C(i), for each gluon 
configuration it is sufhcient to determine the set of quark propagators from 
an arbitrary lattice point to the origin, {S'(a;,0)}. This is achieved by using 
efficient algorithms for the inversion of sparse matrices. 1 now assume that 
these sets of propagators have been generated for many gluon configurations 
and that C{t) has been evaluated. 




(25) 



X 




(26) 
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Figure 1: _Diagramatic representation of the two-point correlation function of the pion, 
eq.(E6|). The lines represent quark propagators in the gluon background field. 

Fig. H shows the pion's correlation function obtained on a 24^ x 48 lattice at 
(3 — 6.2 (for which the inverse lattice spacing is about 2.9 GeV), using 60 gauge 
field configurations (and the SW fermion action, which is an "improvement" 
on the Wilson action as explained in section . The light quarks are a little 
lighter than the physical mass of strange quark. In the same figure is shown 
the fit to a single state. The fit allows for the contribution of the backward 
propagating meson (which arises because, with periodic boundary conditions, 
t is both greater than and smaller than 0, so that there are two contributions 
for any t), so that the exponential is replaced by a hyperbolic cosine (which 
makes a difference close to the centre of the lattice). 

Following the discussion above we need to establish whether we can isolate 
the ground state, i.e. whether the correlation functon is well represented by a 
single exponential over some interval at larger times. To this end it is useful 
to define the effective mass by 



so that if the correlation function is well represented by a single exponential, 
then nics will be almost independent of the time t. The criteria for the con- 
struction of good lattice interpolating operators, such that the overlap with the 
lightest state is enhanced, and a plateau in the plot of rrics (t) as a function of 
t is reached at relatively small times, is an area under intensive investigation. 
In order to determine /^r however, we need to determine the matrix element of 
a local axial current. The effective mass plot for the 7r-mesons, now for light 
quark masses which are a little heavier than the strange quark, obtained by 
the UKQCD collaboration with local interpolating operators are presented in 



Because of finite volume effects it is not possible to perform the computa- 
tions with physical light-quark masses, rather one must extrapolate the results 
obtained for masses around that of the strange quark to the chiral (m^ = 0) 



mcs{t) = In 



(27) 



fig- I 
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Figure 2: Correlation function of the light pscudoscalar meson { "pioir^ with quarks which 
are a little lighter than the strange quark (from ref. l3). 

limit. This requires the knowledge of the value of the K-parameter at the chiral 
limit (kc)- To determine Kc, we use PCAC and plot as a function of 1/k 
(and hence as a linear function of the quark mass), and determine the value of 
K at which = (or takes its physical value) . In fig. ^(a) we show such a 
plot, again from the UKQCD collaboration's simulation at /? = 6.2. In fact it 
can be shown using chiral perturbation theory that the PCAC relation stating 

Sat TO^ is proportional to the quark mass is violated in quenched simulations 
liJ. However this effect only occurs at very low quark masses which is difficult 
to simulate without encountering finite volume effects. Recently Kim and Sin- 
clair have demonstrated this violation of the PCAC in a numerical simulation 
on very large lattices E-3. 

Having determined Kc we now extrapolate the results for /^r to the chiral 
limit. Generally, in the absence of any theoretical guidance, one assumes a 
linear dependence on the mass of the light quark in this extrapolation. Now 
is a dimensionful quantity, which we determine in lattice units, i.e. we compute 
a/jr, where a is the lattice spacing. Thus to make a prediction for we need 
to know the value of a. In lattice computatons this is done implicitly using 
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Figure 3: The effective mass plot of the light pscudoscalar meson ("pioirt with quarks 
which now are a little heavier than the strange quark (from ref.113). 



dimensional transmutation. Instead of fixing a, and then using some physical 
quantity to determine the bare coupling constant go{a) (which is a genuine 
parameter of QCD), one fixes go{a) (or equivalently f3) and determines a. For 
example one may use the mass of the p-meson to determine a (an extrapolation 
to the chiral limit of the UKQCD results, at /? = 6.2, is presented in fig. |^(b)). 
My summary of the values of the lattice spacings for two commonly used values 
of P (obtained by using light quark quantities to set the scale) is 

a-i(/? = 6.0) = 2.0(2) GeV (28) 
a-\(3^ 6.2) = 2.9(3) GeV. (29) 

The errors are largely due to the variation in the results obtained using different 
physical quantities to set the scale, and hence are at least partially due to 
quenching. An inverse lattice spacing of 2 GeV corresponds to a spacing of 
0.1 fm. 

In this way one obtains the matrix element of the axial current (or similarly 
other matrix elements) from lattice simulations. The matrix element is of a 
bare operator which has been defined in QCD with a lattice regularization, 
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3.49 3.50 3.51 3.52 3.53 3.54 
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Figure 4: The values of (a) the pion mass squared, and (b) thM-ho mass plotted as a linear 
function of the quark mass (from ref. t3). 



and a specific lattice action. The ultraviolet properties of such an operator 
are of course different from those defined in the continuum. However, since 
it is the short-distance behaviour which is different, the relation between bare 
lattice operators defined with an inverse lattice spacing 3> Aqcd and 
renormalized continuum operators defined at a scale fi 3> Aqcd , can be derived 
using perturbation theory ||. 

As an example of the results obtained for /tt , the UKQCD collaboration, 
from the simulation which was used to illustrate the calculational steps above, 
quote fTr/rUp — 0.14(1), to be compared to the experimental result of 0.17. 

2 Lecture 2 Simulations in The Heavy Quark Effective Theory 

In this lecture I will apply the formalism reviewed above to lattice simulations 
in the HQET, using the evaluation of the decay constant of the i?-meson (fs) 
as a case study. Matthias Neubert has, at this school, already descibed the 
formalism jXif the Heavy Quark Effective Theory (HQET) and explained its 
usefulness In this framework, predictions for physical quantities in heavy 
quark physics are presented in terms of a series in inverse powers of the mass 

For some quantities this matching between the lattice and continuum operators rep- 
resents the largest source of theoretical upoer*ainty. Non-perturbative methods to evaluate 
the matching factors are being developed EjilZI. 
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of the heavy quark. The predictions are given in terms of matrix elements 
of operators containing the field of the heavy quark, and lattice simulations 
provide the opportunity for the evaluation of these matrix elements. 

Before starting the presentation however, I would like to discuss why it 
useful to perform lattice computations in the HQET, rather than to evaluate 
quantities such as Jb directly in QCD. The reason is that the sizes of the 
lattices and lattice spacings which we can be used are limited by the avail- 
able computing resources, and current quenched simulations are performed 
with inverse lattice spacings typically in the range of 2-3.5 GeV. Thus the 
lattice spacings are larger than the Compton wavelength of the 6-quark, and 
we cannot meaningfully study its propagation directly in QCD. Since in lattice 
computations we have the luxury of varying the masses of the quarks, one ap- 
proach to the study of i3-physics on the lattice is to perform the computations 
with several smaller values of the heavy quark mass (typically in the range of 
that of the charm quark mass), and to extrapolate the results to the &-quark[[ 
Of course such an extrapolation introduces uncertainties, and much effort is 
being spent on reducing the lattice artefacts (errors due to the finite size of the 
lattice spacing) by introducing improved actions and operators, in order to be 
able to perform the simulations with larger masses. 

Simulations in the HQET, provide the results in the infinite mass limit, 
and in principle, also allow for a systematic program of evaluating corrections 
to this limit. The corresponding results can be used as a check of the extrap- 
olations of the results obtained with propagating heavy quarks, as described 
above. Since the mass of the heavy quark does not enter directly into simu- 
lations in the HQET (having been scaled out), lattice artefacts are, a priori, 
expected to be as small as for quantities in light quark physics. 



2.1 The Heavy Quark Effective Theory on the Lattice 



The Lagrangian for the HQET, at zero three-velocity (i.e. in the static limit) 
in Euclidean space takes the form '|: 



C = hD^—-^h 



(30) 



If we define the lattice covariant derivative by 

Di f(x, t) = - \f{x, t) - Ulix, t - a)f{x, t - a) 
a L 



(31) 



"^In the following I shall refer to such simulations as ones with "propagating" heavy quarks, 
in distinction to the "static" quarks of the HQET. 

''The construction of the HQET irt-^jtjfllidean space, with -C ^ is subtle, requiring an 
interpretation of S^''^\x — vt), see refs.li2rE3 
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Figure 5: Quark flow diagram representing the decay of a B-meson into leptons. 



then the propagator of the heavy quark h in the background field configuration 
{U} is given by 

Sl"\x,h(),0) ^ d'-^\x)9{t)ul{6,t - a)Ul{6,t - 2a) ■ ■ ■ Ul{6,0) , (32) 

with S"^*^^ (0, 0; 0, 0) = 1. Other choices of the lattice covariant derivative 
correspond to propagators with diffierent boundary conditions at the origin. 
From eq. ( |32| ) we see that the propagator of the heavy quark is given directly in 
terms of the link variables, and so no inversion of a sparse matrix is necessary 



(as is the case for the quark propagator in QCD, see subsection 1.3 above) 



2.2 Calculation of fs in the Static Limit 

The non-perturbative strong interaction effects in the leptonic decay amplitude 
of the _B-meson are contained in a single parameter, the decay constant fs- The 
process is illustrated in the diagram of fig. ^. To determine t his parameter. 



in analogy to the computation of /tt discussed in subsection 1.3 above, we 
compute the correlation function: 

C2{t) = X^(0|/i(:.)7V9(^)9(0)7VM0)|0) (33) 

X 

= -Y, {Tries' Si{x,0)^S'SHiO,x)]) (34) 



2Mb 

2 f2 

e-^'* + -- - , (36) 



1 \ JbMb ^^£t 



^stat J 2 

where Si and Su are the light and heavy quark propagators respectively. I now 
explain the features of the calculation represented by eqs.(|3^)-(|3^) in some 
more detail: 
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What is the interpretation of £?: E itself is not a physical quantity, 
indeed it diverges linearly with the inverse lattice spacing. At the tree level, 
it is just the difference of the mass of the pseudoscalar meson and that of 
the quark. When higher order corrections are included the interpretation of 
£ becomes more subtle, and will be explained in the third lecture below. In 
the lattice theory, with the hard ultraviolet cut-ofF a~^, quantum corrections 
generate a mass term in the effective theory, even if the bare action does not 
have such a term (as in eq.^O|) ). This mass term is of 0{as/a). For now I will 
just note that from the measured value oi £, it is possible to derive the value 
of any short-distance mass of the heavy quark (such as the MS mass). 

Heavy Quark Scaling Law: Since the action of the HQET is independent 
of the mass of the heavy quark, so is the correlation function (|3^). Thus, from 
eq.(|3^) we can deduce the well known scaling law that the decay constant of 
a heavy pseudoscalar meson decreases like the square-root of its mass. This 
scaling law is modified by logarithmic corrections contained in the matching 
factor Z^'*** as explained below. 

"Smearing": In practice it is jdiffipult to isolate the contribution of the 
ground state in these calculations EilE3 using local interpolating operators. To 
enhance the relative contribution of the lightest state, various extended (or 
"smeared") interpolating operators are used. For example one might smear 
the field in the axial current by 

/i^(x,i) = ^/(2',y)My,i) (37) 

y 

in the Coulomb gauge say. Different choices of the smearing function / are 
used by the various groups carrying out such simulations. The Fermilab group 
in particular have been stressing the possibility of learning about the wave 
function oLthe heavy meson by studying its overlap with different extended 
operatorsLj. However, it should be remembered that Jb is obtained from the 
matrix element of the local (i.e. not smeared) axial current. To determine this, 
one has to compute the correlation functions of two smeared currents (from 
which one obtains the matrix element of the smeared current as in cq.(^5|)) 
and of a smeared and a local current (from which one is then able to deduce 
the matrix element of the local current). A strong check that the ground state 
has been isolated is provided by verifying that the matrix element of the local 
current is independent of the smearing function / used in the calculation 
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Matching: The quantity which one obtains directly in lattice simulations is 
the matrix element of the bare axial current in the HQET, defined with the 
lattice regularization. In order to obtain /b one needs the matrix element of 
the axial current in QCD. The difference between the two is an ultraviolet 
effect and can be calculated in perturbation theory. Although not necessay, 
this calculation is usually performed in two steps 

^latt, HQET(^) ^ ^MS, HQET(^) _^ ^QCD (33) 

In the first of these steps lattice perturbation theory is used to obtain the 
matrix element of the axial current in the HQET defined in some continuum 
renormalization scheme, such as the MS scheme, from the matrix element of the 
bare current in the lattice theory. In the second step the matrix element of the 
physical (QCD) current (and hence fs) is obtained from that in the effective 
theory. Z^*^* in eq.(|3^) is the combined matching factor. It depends loga- 
rithmically on the lattice spacing, being proportional to as{a)'--~^^^°\ where 
Po = 11 — 2/3Nf is the coefficient of the first term in the /^-function, and Nf 
is the number of light-quark flavours. 

In these lectures it is not possible for me to describe the technical details 
of lattice perturbation theory. The lattice action for the HQET in eq.(pO|) 
and the QCD action in eq. (pT[) , with the derivatives defined as differences, 
lead to Feynman rules where the momentum dependence is given in terms 
of trigonometric functions. For example the quark propagator for the Wilson 
action is given by the expression in eq. ( p^ ) . This makes the analytic evaluation 
of Feynman diagrams in lattice perturbation theory prohibitatively difficult in 
general, and instead the integrals are evaluated numerically. 

Lepage and Mackenzie have argued convincingly that the bare lattice cou- 
pling 59(0) is a poor expansion parameter for lattice perturbation theory, i.e. 
that withjthis choice the higher order corrections have unnecessarily large co- 
efficientso. As an example, consider the one-loop tadpole correction to the 
plaquette variable, illustrated in fig. |^. Numerically it turns out that such 
tadpole contributions, which are common to many quantities, are large. In 
particular they will appear each time that there is a link variable (from which 
the tadpole can emanate). Let 

ut^^{Ti-V,,.) (39) 

where the plaquette V has been defined in eq. (|l2|) . Thus one might expect 
that at least part of the large higher order corrections are given by uq for each 
link variable present in the quantity being studied. If this is the case then we 
can replace each link variable by U^/uq, compensating for this by different 
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Figure 6: Representation of a one loop tadpole diagram contributing to the plaquette vari- 
able. The diagram is generated by expanding one of the link variables to 0(gg). 



perturbative matching factors. For the gluon action one would rewrite V / 
7'/(g^Ug), with — g^/uQ. It might therefore be expected that is a better 
expansion parameter, and this appear to be the case (for numerical evideace of 
this feature, and of a general discussion of "tadpole resummation" see ref.cJ) . It 
would, nevetheless, be very welcome to have a non-pertubative determination 
of the matching factor. For simulations with the Wilson action for the light 
quark, Z'^ ~ 0.7 for typical beta values used in current simulations, and for the 
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"improved" SW-action described in subsection |2.4| , ~ 0.8. These values 
are obtained from one-loop perturbation theory. 



2.3 Results 

Before presenting some results obtained from simulations in the heavy quark 
effective theory, let us recall precisely what it is that we are calculating. The 
decay constant (/p) of a heavy pseudoscalar meson {P) satisfies the scaling 
law: 

/p = -^|a;2/^nMp)(l + 0(a,)) +0(^)| , (40) 

where Mp is the mass of the meson. The factor of as{Mp)~^^^° arises from 
the matching of the HQET with QCD, and in our case above (for the _B-meson) 
it is included in Z^^^. From lattice simulations in the HQET, we obtain A, 
the normalisation of the leading term. In order to compute the 0{1/Mp) 
corrections, it is necessary to compute the matrix elements of the dimension-4 
operators which arise in the expansion of the axial current in QCD in terms of 
operators of the HQET. It is also necessary to include the 0{l/Mp) corrections 
to the heavy auark action (|30|). Such calculations are only just beginning to 
be performed c3. I denote by f^^^, the value of fs obtained by dropping the 
1/Mp corrections. The UKQCD collaboration obtain 

/l*'^* = 266 1 2o(stat) i 27(syst) MeV , (41) 



17 



where the quoted systematic error is due to the uncertainty in the value of the 
lattice spacing (which enters as c?^^ in this calculation). For a recent review 
of the status of computations of /^*^* see ref. □. I will delay the interpreta- 
tion of this result until after a discussion of the computations performed with 



propagating heavy quarks in subsection 2.4. 

The calculation of Z^'*^' is performed in a theory without any explicit large 
masses. A priori, one would therefore expect that the systematic errors should 
be similar to those in light quark physics. By repeating this calculation at sev- 
eral values, of /3, the Fermilab group have claimed that there are large lattice 
artefactso. This conclusion relies heavily on the results at /3 = 5.7 however, 
and a detailed study of the data of the APE collaboration leads to the copflu- 
sion that the results for are independent of /? (i.e. of a) for /? > 6.0E3. 

2.^ Calculations of fs performed with Propagating Heavy Quarks 

In this subsection I return to simulations performed in QCD. Now, as explained 
above, the 6-quark is too heavy to use, so the calculations of decay constants 
and other physical quantities are performed with quark masses in the region 
of that of the charmed quark. Even so, the largest source of systematic uncer- 
tainty are discretisation errors which are of O^mga) for a heavy quark with 
mass mq. For example, at /3 = 6.2 the bare mass of the charm quark, rrichfij 
is about 0.35 in lattice units, where 

1 / 1 1 \ 

mch,oa = - I I (42) 

and Kch is the K-parameter corresponding to a charm quark. 

"Improvement": In order to reduce these discretisation errors much ef- 
fort is being invested to develop the systematic "improvement" technique of 
SymanzikllJ. This technique involves the modification of the lattice action and 
operators in such a way that the_discretisation errors are formally reduced. As 
an example consider the actioncS: 

^QCD ~ '5'gluon + 'S'^^ AS* (43) 

where S'giuon and S^r are the Wilson gluon and quark actions as defined in 



subsection 1.1 and 



AS"" = -2Ka^ V ■^\ip{x)Uf,{x)U^{x + afi)^p{x + 2afl)+ 

^ — ^ 8(2 

^pix + 2afi)Ul{x + ajl)Ul{x)ij{x) - 2{l){x)il){x)^ . (44) 
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Now expanding Sqqj^ as a series in the lattice spacing we find: 

Sqcd = S'^cD + 0(«') , (45) 

where the superscript "cont" represents "continuum" . A similar expansion for 
the Wilson action, eg. ([l8|) , gives discretisation errors of 0{a) rather than of 
0{a?). When radiati|Vfi corrections are included one finds that the discretisation 
errors are of 0{asa)^, and so the use of the action ( ^3|) has formally reduced 
the lattice artefacts from 0{a) to 0{asa). One can also readily check that the 
problem of fermion doubling has not been reintroduced by the addition of the 
term AS*^^. 

The action in eq.(|4^) contains next-to- nearest neighbour interactions, in 
addition to the nearest-neighbour interactions of the Wilson theory. By chang- 
ing (the fermionic) variables in the functional integral, it is also possible to 
constmct an equivalent improved theory with nearest neighbour interactions 
onlyE3 

+ S^~ AS' (46) 

with 

= ^go E ^i^W^F^M^) . (47) 

This lattice action, which was first proposed by Sheihkoleslami and Wohlert, 
is called the SW-action. In eq. (^7|), F^^ is a lattice definition of the field 
strength tensor. Several groups performing large scale numerical simulations 
are using this action, and the results from the UKQCD collaboration which I 
am using to illustrate this talk were obtained with this formulation of QCD. 
Even with this action however, in simulations containing the charm quark, one 
might expect 10% errors due to finite-a effects. 

It is also possible to try to construct an action by subtracting AS' with a 
coeflacient which is different from unity in order to reduceJihe lattice artefacts 
still further. This is sometimes done at the one-looD level Ell or by using mean 
field theory estimates or non-perturbatively Ell. Eiiially let me mention 
the exciting program of constructing "Perfect Actions" a, in which, ideally, by 
using the renormalization group transformations on a fixed-point action, all the 
discretisation errors would be removed (although in practice one may have to 
settle for a substantial reduction of the artefacts, rather than for "perfection"). 

Results: My summary of the results obtained with propagating quarks is: 

fo = 200±30MeV (48) 
fs = 180±40MeV, (49) 
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where the result for the _B-meson was obtamed, of course, by extrapolation 
to the 6-quark of results obtained with lighter quarks. This extrapolation 
is illustrated in fig. 0, where the results from the UKQCD collaboration are 
presented for the scaling quantity 

$(Afp) = fp {as{Mp)/as{MB)f'^'' , (50) 

as a function of l/AfpEl, where all quantities are given in lattice units |. Here 
P represents a heavy-light pseudoscalar meson. There are twelve "measured" 
points, corresponding to three light quark masses (represented by the open 
triangle (heaviest), circle and square(lightest)) and four heavy quark masses 
(represented, for each light quark mass, by the same open symbol). The full 
points represent the values obtained by extrapolation of the light-quark masses 
to the chiral limit, for each heavy quark mass. If there were no \/Mp and 
higher corrections then, for each choice of the light quark mass, the points 
would lie on a horizontal line. Clearly there are significant, negative, 0{1/Mp) 
corrections. The violation of the heavy quark scaling law is of the correct sign 
to match up with the static result, although UKQCD find a (small) discrepancy 
between their result for and the value obtained by extrapolating the 

results computed at finite values of Mp. Given the different systematic errors 
in the computations with static and with propagating heavy quarks, the near 
consistency of the results is very pleasing. I would summarize the best estimate 
of Fb from lattice simulations as being 190(40) MeV. 



Test of the Heavy Quark Symmetry: I would like to mention one im- 
portant test of the lattice calculations withjjropagating heavy quarks, which 
is provided by the heavy quark symmetry E3. In the infinite mass limit the 
decay constants of the pseudoscalar {fp) and vector {fv) mesons are related 



as M ^ oo, where fv is defined by 

{Q\V,mV) = e,-^ . (52) 
Jv 

In eq.(^), is the vector current with the appropriate flavour quantum num- 
bers, and is the polarisation vector of the vector meson V . The perturbative 

"^The simulation was performed at /3 = 6.2 for which the inverse lattice spacing is about 
2.9 GeV. 



20 




Figure 7: The behaviour of the scaling quantity <I> with l/Mp. The solid line represents the 
hnear fit to the chirally extrapolated points using the three heaviest meson masses, whereas 
the dashed curve results from a quadratic fit to all four. 



matching factor in eq.(|5l|), has been calculated to one- loop order. In fig. |5| I 
present the results of the ratio in eq.(|5|) as a function of 1/M, where for M 
I have taken the spin averaged meson mass {M — (3 My + Afp)/4) M. It 
can be seen that the results obtained at finite values of the mass of the heavy 
quark differ substantially from the asymptotic value, but the extrapolation to 
M = oo is in good agreement. 

2.5 Other Physical Quantities in B -Physics 

To conclude this lecture I will briefiy mention some other quantities in heavy 
quark physics which have been, or are being, computed in simulations in the 
HQET. 

fBa/fBd- -By studying the dependence on the results for /b on the mass of the 
light quark, one obtains a good estimate of the ratio of the decay constants 
of the Bs and Bd-mesons. With both static and propagating heavy quarks 
one typically finds that fs, is about 10-20% larger than /s^. This can be 
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Figure 8: The quantity U{M) plotted against the inverse spin-averaged mass. Linear and 
quadratic fits are represented by the solid and dashed curves respectively. Also shown are 
the statistical errors of the extrapolation to the infinite mass limit. 

understood as being due to the fact that the Bs meson is more compact than 
Bd, with a larger wave- function at the origin. 

Heavy Baryons: The detailed study of heavy baryons is just beginning in 
lattice simulations. In the HQET it is possible to evaluate mass differences, 
such as that between the Ai, baryon and i?-meson. The UKQCD collaboration 
find for this quantity EZl: 

Ma, - Mb = 420 t'^^o ^ 30 MeV , (53) 

to be compared with the experimental result of 346 MeV with an erroiLX)f 
about 7 MeV (this is based on the new result from the CDF collaboration 
It must be remembered that the result in eq.(p3|) represents the value in the 
infinite mass limit. A comprehensive study of the spectrum of heavy baryons 
using propagating heavy quarks, can be found in ref.E^, where the Ab-B mass 
splitting was found to be 359 ± 50 ± 27 MeV. 
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Bb: One of the most important quantities in the phenomenology of _B-physics 
is the parameter Bb, which together with fs contains the strong interaction 
effects in B° ~ B mixing. It is defined by 

where Ol is the AB — 2 operator 

Ol(a^) = - l'')qbr{l - l')q ■ (55) 

The first factor on the right hand side of eq.(p4|) is introduced to cancel the 
H (i.e. renormalization scale) dependence of the matrix element at leading 
logarithmic order. 

In order to determine the physical amplitude for the A_B — 2 transition, it 
is necessary to write the corresponding QCD operator, renormalized in some 
continuum scheme, in terms of the bare operators of the HQET defined us- 
ing the lattice regularization. Because of the breaking of chiral symmetry in 
theories based on the Wilson formulation of lattice fermions, this requires the 
evaluation of the matrix elements of 3 other A_B — 2 operators in addition to 

Or - (&7p(1+7')9)(^7m(1+7')'z) (56) 
On = (&7p(l-7')'?) (^7m(1+7')'z) + 

2 (6(1 - 75)q) (5(1 + 7')(Z) + 2 (5(1 + 75)g) (6(1 - ^^)q) (57) 

Os = (6(1 - 7'')'?) (5(1 - 7')?) , (58) 

where the QCD operator O^ is given by 

0?^^(mfc) = Zl Oi"" + Zr O'^'' + Zn OJ^" + Zs O^"" . (59) 

The matching coefficients Zi, {i — L, R, N, S}, have been calculated at one 
loop level for the Wilson and SW actions in references e3 andO respectively. 
The one-loop contribution to Zl is large, for example at /3 = 6.2 with the SW 
action (corresponding to the UKQCD simulation), Zl — 0.55. This underlines 
further the necessity for a non-perturbative determination of the matching 
coefficients. 

The matrix elements of the lattice four-quark operators, normalised by 
the factor 8/3/^m^ as in eg. (^^ , can readily be obtained from the ratio of 
correlation functions: 

nSS(. . N _ Kf^{ti,t2) 
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Figure 9: Diagramatic representation of the ratio defined in eq.(^). The double (single) 
lines represent the propagators of the heavy (light) quark. The labels S and L indicate 
whether smeared or local operators have been inserted at the corresponding point. 



where K^^ , with i = L,R, N, S, is the three point function 

KfHti,t2) ^ ^ (0|{At'(fi,-ti)Or(0)At'(f2,i2)} |0)(61) 

Xi,X2 

C^^ is the two-point correlation function of two axial currents, 

X 

and the superscripts L and S stand for local and smeared respectively. The 
smeared currents are introduced to enhance the overlap with the ground state, 
and the local ones because it is from these that we are able to obtain the decay 
constant fs- The ratio Ri is represented by the diagram in fig.|^. 

The UKQCp, collaboration finds the following results from their simula- 
tions at /3 = 6.2E3: 

Bb, = 1.02 t ^(stat) 1 ^(syst) (64) 
Bb, = 1.04 t 5(stat) 1 ?(syst) (65) 
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A discussiop_.of the phenomenological consequences of these results can be 
found in refO. A preliminary result which is somewhat higheiL.(by about 25%) 
than that in eqs.(|6^ and ( |65| ) was recently presented in ref.LJ. 

Relatively few calculatiDnsiiave been performed for the B-parameter with 
propagating heavy quarksCJuCj, the most recent of which quotes 

ai^/'^°)(5 GeV) Bb = 0.90(5) and 0.84(6) (66) 

obtained using the Wilson ferniion action at /3 = 6.1 and 6.3 respectively^ Q 
Thus this result obtained with propagating Wilson quarks is a little larger 
(about 20% larger) than that obtained with static quarks in eq. (|64[) . It is 
too early to be sure that this difference is a physical effect, rather than an 
artefact of lattice systematics. In particular it will be interestixLflj_tp compare 



the calculations of the matching factors in the recent studipaoo when t 



details of these calculations are published, with those of ref.tZfE^. In referenceE 
there is some evidence that the 0(1/toq) corrections to the result obtained in 
the static limit may be negative, so that the _B-parameter increases slowly as 
the mass of the heavy quark increases, which is the opposite conclusion one 
would draw from comparing the results in eqs.(|6^ and (|6^), with that of the 
JLQCD collaboration with propagating quarks. 



The Hyperfine Splitting, Mb- — Mb- As a final example of quantities 
which can be evaluated in lattice simulations, consider the mass difference 
between the vector and pseudoscalar heavy-light mesons. In the static limit 
these mesons are degenerate, and their mass difference is of ©(l/mg), where 
toq is the mass of the heavy quark. The hyperfine splitting is given by the 
matrix element of the chromomagnetic operator which appears at 0{l/mQ) 
in the action. The results of lattice calculations of the hyperfine splitting 
give values which axe, cprisiderably lower (by factor of 2 or so) than the 
experimental result c^EZlca. For example in ref.1221 it was found that 

M|. - Ml = 0.281 ; \l (stat) t 3? (syst) GeV^ , (67) 

where the systematic error represents the uncertainty in the value of the lattice 
spacing only. The experimental result for this quantity is 0.488(6) GeV^ell. 

It has been observed for a considerable time that the values of the hyperfine 
splitting in heavy mesons (including heavy quarkonia) are found to be too small 
in lattice simulations with propagating heavy quarks. It is understood that 
this physical quantity has particularly large discretisation errors, and indeed 

f Thus the JLQCD group quote their results in terms of the matrix element *anormalized 
at 5 GeV in the MS scheme, and not the renormalization group improved oneCj. 
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as one "improves" the action, the disC|Lepency between the lattice resuhs and 
experimental measurement decreases C3t3. It is peculiar however that the 
result obtained with the HQET is so different from the experimental value. 
In the HQET the discretisation errors are all of ©(aAqco), since there is no 
dependence on the mass of the heavy quark. Of course quenching may be 
partly responsible for the discrepency, and in addition the matching factor 
(between the lattice operator in the HQET and that in QCD) is found to have 
very large one-loop corrections (the factor, which is 1 at the tree level, is found 
to be about 1.8 at the one- loop level). This is another example of the necessity 
of computing these matching factors non-perturbatively. 



3 Lecture 3 — Renormalons and Power Divergences in Lattice Sim- 
ulations 



In the final lecture I will briefly review the status of a subject which is currently 
very topical, not just in lattice simulations, but in phencjmenological studies 
of QCD in general. This is the problem of renormalons e2I, which arises when 
one is interested in computing power corrections to hard scattering and decay 
processes. For example, in the discussion above this problem would arise if one 
were to compute the 0(1/toq) corrections to the decay constant of a heavy 
meson. Other important examples include the evaluation of higher twist effects 
in deep- inelastic structure functions and the contributions of the giuon conden- 
sate to quantities which are measured in e+-e^ annihilation. The discussion in 
this lecture will be restricted to those processes for which an operator product 
expansion exists. 



3. 1 Overview 

I will now present the principal results concerning the appearance and cancel- 
lation of renormalon singularities in operator product expansions, using the 
Heavy Quark Expansion as an illustration. Consider the evaluation of the ma- 
trix element of some local QCD operator containing one or more heavy quark 
fields. We call this operator O^'-'^. For example if we wish to evaluate the 



leptonic decay constant Jb as in section 2.2, then O^*^^ would be the axial cur- 
rent b-f^j^q, where q represents the field of the light quark. Using the HQET 
we expand O^'-'^ in inverse powers of the mass of the heavy quark (?tiq): 



OQCD^Cr(mQ/M)OfQ^^(^) + 
1 



— C2(™Q//i) Of'^^'^i,,) + Oil/ml) (68) 
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where /i is the scale at which the operators of the HQET, o^Q^""" and O^^^"^ , 
are renormaUzed. We consider here the simple situation for which there is 
a single operator in each of the first two orders of the expansion, but the 
discussion below can be easily generalized to the case in which there are more 
operators. 

When evaluating matrix elements of O^*^^ beyond the leading order in the 
1 /mq expansion, in addition to the corresponding higher dimension operators 
on the r.h.s. of eg. (|68|) , it is also necessary to keep higher order terms in the 
heavy quark action, i.e. by replacing the Dirac term in the QCD action by 

Q{ij ■ D - mQ)Q — > hy{iv ■ D)hy + hy{iDfhy 

+ ^^hyC7^0F"^hy + 0{l/ml) , (69) 

ZlTlQ Z 

where Q and hy represent the fields of the heavy quark in QCD and the HQET 
respectively, and v is the quark's four-velocity. In the discussion below we will 
take the quark to be at rest, and will denote the corresponding heavy quark 
field by h (without any subscript). Cmag is a constant determined by matching 
the effective theory onto the full one (QCD). The corresponding constanl^for 
the kinetic term {h{iD)^h) is equal to one by reparametrization invarianceEata. 
Throughout the discussion below it is implied that the action of the HQET 
contains sufficiently many terms for the precision of the calculation. 

In general the QCD operator OQ'^^ also requires renormalization and is 
defined at some scale M. Unless specifically needed I will suppress the depen- 
dence on M in O^^'^ and in the coefficient functions C^. 

I will now summarize the main points w hich T wish to make in this lecture, 
a more detailed account can be found in ref. Hd'Ej: 

• If we restrict the calculation to the leading order in I/toq, i.e. if we 
neglect the ©(l/mg) terms, then the perturbation series for the coef- 
ficient function Ci diverges, and is not Borel summable. This is due 
to a singularity in the Borel transform of the perturbation series, called 
an infra-red renormalon. This implies that there is an ambiguity in the 
evaluation of Ci, coming from the different possible ways of defining the 
sum of the series. This ambiguity is of ©(l/mg). 

• Formally therefore, we should not include the 0(1/toq) corrections in 
eq. ( |6^ ) until we have computed sufficiently many terms in the perturba- 
tion series for Ci to control its divergent behaviour. 

• For this talk I will restrict the discussion to problems for which the matrix 
elements of 0^'~^^ have no renormalon ambiguities. In general such non- 
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perturbative effects do exist, and appear on the right hand side of eq.(p8[) 
in the matrix elements of the operators of the HQET. They do not affect 
the coefficient functions. 

In renormahzation schemes based on the dimensional regularization of 
ultraviolet divergences in the HQET, such as the MS renormahzation 
scheme, the matrix elements of O^^^""" are also not Borel summable 
in perturbation theory, due to an ultraviolet renormalon singularity in 
their Borel transforms. The corresponding ambiguities in the matrix 
elements of O^^^^ cancel those in Ci. Predictions for physical quan- 
tities cannot have any ambiguity, of course, but the matrix elements of 
higher-dimensional, or higher-twist, operators are, in general, not defined 
uniquely by the MS procedure. 

If a hard ultraviolet cut-off is used, such as the lattice spacing in the 
lattice HQET, then the matrix elements of O^'^^^ do not have am- 
biguities due to ultraviolet renormalons. Indeed the matrix elements 
can be computed (unambiguously) in lattice simulations, in contrast to 
those in the MS scheme, which cannot be computed directly using some 
non-perturbative technique. This, in turn, implies that the coefficient 
functions Ci do not contain ambiguities due to infra-red renormalons. 
In the lattice theory it is natural to present the discussion in terms of 
bare operators in the effective theory, defined with the lattice spacing a 
as the ultraviolet cut-off. I will assume here that mga ^ 1. Of course if 
the inverse lattice spacing was much smaller than the heavy quark mass 
then there would be no need to use the HQET. 

The absence of renormalon ambiguities in Ci in the lattice theory is 
due to a cancellation between terms which, in any order of perturba- 
tion theory, are of different order in the l/mg expansion, i.e. between 
terms which behave logarithmically with mga and non- leading ones of 
0{l/mQa). 

In the lattice theory, matrix elements of higher dimensional operators 
(such as O^'^^^) diverge as inverse powers of the lattice spacing, and are 
hence manifestly unphysical. For example, for the matrix element of the 
kinetic energy operator, 

{H\hD^h\H) ^ Oil/a"^) (70) 

where H represents the heavy hadron. This can readily be confirmed in 
perturbation theory. 
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• The subtraction of the power divergences in perturbation theory intro- 
duces renormalon ambiguities. Thus for example, the perturbation series 
of the terms which diverge quadraticaUy, i.e. those which are 0{\/a?), 
in eq.([70|) contains a renormalon ambiguity. 

• The pprtiyhatinn scrics for the pole mass has a renormalon ambiguity of 
0(AQCD)oEa. This is not the case for short-distance definitions of the 
heavy quark mass, such as 

ruQ ^ mf{m^) . (71) 

• It is possible to compute mq (and other short-distance definitions of 
the mass), using only simulations in the HQET. The quantity which is 
computed directly in lattice simulations is the bare binding energy, £ 
in eq. (|35|) . E is also not a physical quantity, diverging linearly with the 
inverse lattice spacing. It's significance can be deduced by matching 
QCD and the lattice formulation of the HQET, from which one finds 

£ — Mp — TTipoic + 5m , (72) 

where Mp is the mass of the pseudoscalar meson ^, mpoic is the pole mass 
of the heavy quark, and 5m the perturbative series of linearly divergent 
terms contributing to. the mass renormalization in the effective theory 
with the action (pO|)E3. The perturbation series for mpoic in terms of m 
has a renormalon singularity, which is cancelled by that in 5m. The linear 
divergence in 5m cancels that in £. Thus using eq.([72[), the measured 
value of S and perturbation theory, one can obtain to. In refE3 it was 
found in this way that 

nib = 4.17 ± 0.05 ± 0.03 GeV + 0{l/mb) . (73) 



• The presence of renormalons in physical quantities for which there is no 
Operator Product Expansion, such as the Drell-Yan process or in event 
shape vadahles in jet physics, is a subject currently under intensive inves- 
tigationEJ'ES. It is hoped that these studies will provide important phe- 
nomenological information about the sub-asymptotic (non-leading twist) 
behaviour of physical quantities. In lattice QCD there are analogous 
questions about the presence of ambiguities in perturbation series for 
quantities which contain power divergences, e.g. is there an ambiguity 
in the perturbative evaluation of the critical mass when using Wilson 
fermions? 

^ Other hadronic states can also be used equally well to determine m. 
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Final Remarks 



I hope that in these lectures I have managed to convey the exciting possibilities 
which lattice simulations offer for the computation of non-perturbative strong 
interaction effects in physical quantities. The evaluation of these quantities 
is important to progress in particle physics. The applications to _B-physics 
which were discussed in the secong lecture, are only a small subset of the 
phenomelogically relevant quajotities which are being computed using "on the 
lattice" (see for example refs.ElBB). At the same time I have tried to explain 
the difhculties which remain, of which a reliable formalism and algorithms for 
the inclusion of light quark loops is, mot probably, the most significant. 
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